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We study the dynamics of a dimer moving on a periodic one-dimensional substrate as a function 
of the initial kinetic energy at zero temperature. The aim is to describe, in a simplified picture, 
the microscopic dynamics of diatomic molecules on periodic surfaces, which is of importance for 
thin film formation and crystal growth. We find a complex behaviour, characterized by a variety of 
dynamical regimes, namely oscillatory, "quasi-diffusive" (chaotic) and drift motion. Parametrically 
resonant excitations of internal vibrations can be induced both by oscillatory and drift motion of the 
centre of mass. For weakly bound dimers a chaotic regime is found for a whole range of velocities 
between two non-chaotic phases at low and high kinetic energy. The chaotic features have been 
monitored by studying the Lyapunov exponents and the power spectra. Moreover, for a short-range 
interaction, the dimer can dissociate due to the parametric excitation of the internal motion. 
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I. INTRODUCTION 



Thin-film growth is a topic of great importance both 
from the theoretical and experimental point of view [1] . 
In order to build a microscopic theory of crystal growth, 
it is fundamental to understand isolated-adatom surface 
diffusion, and indeed a great amount of work has been de- 
voted to the monomer case during the last decades [2-6] . 
Once individual atoms are adsorbed on a surface, they 
can meet each other thus forming dimers. In spite of 
this very simple morphology, dimers display a peculiar 
and intriguing diffusive behaviour [7-14]. In this paper 
we show that the strong nonlinearities arising during the 
motion of dimers on a periodic substrate make dimer dy- 
namics a very complex phenomenon. Here, we study the 
problem of a one-dimensional dimer moving on a rigid 
periodic substrate at T = 0. A one-dimensional model 
is relevant since one-dimensional dimer diffusion occurs 
in real systems, in particular along steps and on chan- 
neled (110) metal surfaces [15]. Although some aspects 
of dimer energetics and dynamics have been recently con- 
sidered (see in particular Refs. [8,9,11]), here we focus on 
the nonlinear microscopic dynamics of this model, which 
gives rise to different regimes depending on the initial 
velocity and on the strength of the dimer interaction. 
In particular, chaotic motion with long jumps is found 
between the oscillatory behaviour at low initial kinetic 
energy and the drift motion at high initial kinetic en- 
ergy. In all these regimes, situations can be found where 
the centre of mass (CM) motion drives excitations of the 
dimer vibrations which can lead to dissociation for short 
range interatomic interactions. Such a complexity in the 
dynamical behaviour at T — is relevant to understand 
the thermal diffusion problem. In particular, the chaotic 



behaviour which dominates for weakly bound dimers can 
account for the non-trivial dependence of the diffusion 
constant on the strength of the interatomic interaction 
reported in Ref. [9]. 

Chaotic motion can occur in nonlinear systems with at 
least three variables. Much studied is the case of systems 
characterized by a single spatial coordinate subjected to 
an external drive [16-20]. Besides, a system of inter- 
acting particles can exhibit chaotic motion with "quasi- 
diffusive" features in one dimension [21-23]. This can 
happen even without the presence of an external drive, 
as we will show below. In view of the sinusoidal potential 
and of the phase space dimension, our model bears some 
resemblance to an undriven double pendulum [24-26]. 
However, in our case, the possibility to perform either 
oscillatory or drift motion leads to the appearance and 
subsequent disappearance of chaos for increasing initial 
kinetic energy. 

In Sect. II we describe our model. Sect. Ill is devoted 
to the discussion of the linearized version of the system, 
in which a semi-analytical treatment can be performed. 
Sect. IV analyses the dynamics of the full system, and in 
Sect. V we discuss the chaotic properties of the system. 
Some concluding remarks are presented in Sect. VI. 



II. MODEL 



We consider a dimer moving on a periodic one- 
dimensional substrate at zero temperature without 
damping. The particle-substrate interaction is modelled 
by a periodic function U : 



U(xi,x 2 ) = Uq(2 — cos(fcKi) — cos(kx2)) 



(1) 
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where Xi represents the spatial coordinate of particle i 
(i = 1,2), 2Uq is the potential barrier per particle and 
k = 2ir/a, a being the substrate lattice constant. Most 
of the results presented here have been obtained using a 
harmonic interatomic potential: 



V(X1,X 2 ) = y (X 2 ~ Xl - if 



(2) 



where K is the force constant and I is the spring equi- 
librium length. We also used short range interatomic 
interactions, such as the Lennard- Jones (LJ) potential 
(see Sect. IV B). 

The equations of motion are: 



mil = K(x2 
mi 2 = K(x\ 



xi-l) 
x 2 + I) 



kilo sin(fca;i) 
kilo sm(kx2) 



As initial conditions we choose 

xcm(0) =x Q = a/2 x C m(0) = v 
x r (0) = x r (0)=0. 

In this way we give an initial kinetic energy to the dimer 
at equilibrium (alternatively one could have chosen to 
give an initial potential energy, i.e. xcm(O) = and 
icm(O) 7^ °/2). The minimum kinetic energy for the 
CM to get out of the potential is Vq = 4[/ if v\ = v 2 = vq 
(vi is the initial velocity of particle i). Below this thresh- 
old value (namely v < ^f4^o) Eq. (6a) coincides with 
that of a classical pendulum for which the amplitude as 
a function of the period of oscillation is known in terms 
of an elliptic integral [27], as shown in Fig. 1(a). 



(3) 



where m is the mass of each particle. We rescale the 
variables in the following way: 



kxi, t = t/r, U a = U /E T , 



kl, K = K/ (Ex 



where Et is a reference energy and r = [m/ (ETk 2 )] 1 / 2 . 
In these units the equations of motion become (in the 
following we will omit the tildes for simplicity) 



Xl 
1 2 



K(x 2 



Xl-l) 

x 2 + 1) 



Uo smxi 
U sin x 2 



In the CM and relative coordinates frame we have 

f x cm = -U Q smx C M cos(x r /2 + 1/2) 

\ x r = — 2Kx r — 2Uq cosxcm sin(x r /2 + 1/2) 



(4) 



(5) 




where xcm — {x\ + x 2 )/2 is the CM coordinate and 
x r = X2 — xi — I is the relative coordinate. We concen- 
trate on the commensurate case in which I = a = 27r, i.e. 
the spring natural length is equal to the period of the sub- 
strate potential. In this situation the minimum energy 
configuration does not depend on K (namely x\ — and 
x 2 = a minimize the total potential energy) and more- 
over a linearization around x r = offers the possibility 
to treat the problem in a semi-analytical way. This has 
the advantage to give a closer insight on the dynamical 
features of this system. 

We perform molecular dynamics simulations, integrat- 
ing the equations of motion (4) using a velocity- Vcrlet 
algorithm, with time step A = 10~ 4 . 



III. DYNAMICS OF THE LINEARIZED SYSTEM 

When x r ~ 0, as at the beginning of the motion start- 
ing from equilibrium, we can linearize in x r the equations 
of motion (5): 
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FIG. 1. Relation between amplitude and frequency of os- 
cillation as given by the solution of Eq. 6a (a), between CM 
frequency of oscillation and CM initial velocity (b) and be- 
tween the effective stretching frequency and CM initial veloc- 
ity (c), obtained by numerical calculations for Uo = 0.6 (solid 
lines). The lower borders of instability regions as a function 
of K are plotted in (d) . The vertical dashed lines in (a) and 
in (b) indicate respectively the frequency of small oscillations 
VUo = 0.774 and the velocity corresponding to the thresh- 
old for drift motion y/4Uo = 1.5492. The horizontal dashed 
lines in (c) and (d) correspond respectively to wo e // = V2K 
(K — 0.1), which is reached when vo > y/AUo, and to the 
energy threshold for drift motion y/AUo- 

The maximum amplitude of oscillation of the CM is 
determined by vq. In this case 



C 

x C m ~ x + — sin(wt). 



(7) 



(8) 



xcm ^ U sin x C M 



-2K{1 - ^ cosxcm)xt- 



(6a) 
(6b) 



with lo =ui(vo) (see Fig. 1(b)). This means that 

cos xcm — A + Bcos(2ojt) 

with A and B depending on v via uj. Inserting (8) into 
the equation of motion for x r (6b) we obtain the equation 
of a parametric oscillator: 



'x r = — Wg e //(1 + hcos(2Lot))x r 



(9) 
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where uioeff = yJ2K — AU^ and h = BUq/uXq^j. It is 
worth noting that the stretching frequency of the free 
dirner y/2K becomes \/2K + Uq in the external poten- 
tial (this is true when the CM is fixed at the equilib- 
rium position). But in Eq. (6a) the motion of x r is fur- 
ther affected by the oscillatory behaviour of the CM and 
its natural frequency changes into u>oeff = \J2K — AU§. 
Conversely, when v > VWo the CM performs a drift 
motion, i.e. 

X C M ^ X + < V C M > t, (10) 



The different curves represent ftasa function of u> for 
different values of K. Note that, at fixed K, h increases 
when u> decreases (e.g. when the amplitude of oscilla- 
tion increases and vo < y/4Uo), but when the CM over- 
comes the barrier, h reaches the constant value Uq/{2K). 
Moreover, the range of frequency in which instability is 
observed is larger for smaller values of K. An example of 
parametric resonance is shown in Fig. 3(a): x r oscillates 
and its amplitude increases exponentially. In Fig. 3(b) 
a blow up of the behaviour of x r and of the drive xqm , 
oscillating at the same frequency, is also shown. 



where < • > denotes time averages and < vcm > is the 
drift CM velocity. By defining w as w =< vcm > /2 
the equation for x r remains in the form (9) with woe/ / = 
V2K and h = U /(2K) (i.e. A = and B = 1). 

The CM motion (either oscillatory or drifting), as con- 
sidered in the linearized Eq. (6a), drives paramctrically 
the internal motion of the dimcr. We establish the in- 
stability ranges by monitoring for which values of the 
initial velocity v$ an exponential increase of x r is found. 
The relation between the initial velocity vq and the fre- 
quency uJoeff is shown in Fig. 1(c) for several values of 
v < y/AUo- In order to understand the energy thresh- 
old for the excitation of parametric resonances, we plot 
the lower borders of the instability regions as a function 
of K in Fig. 1(d) when the total initial kinetic energy 
E® in — Vq is less than the potential barrier 4U . In this 
way we can identify a critical value K — K c above which 
the parametric resonance can be excited only if the CM 
performs a drift motion (e.g. E® in > 4C/ )- It turns out 
that K c ~ 0.3 for Uo = 0.6. Since woe// is determined 
by K, by considering different values of K we can con- 
struct the standard picture for parametric instabilities 
relating /i to w. One can recognize the main resonance 
for u) = oJoeff (Fig. 2). The boundaries of the region of 
instability are given by the stars in Fig. 2. 
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FIG. 2. Relation between frequency and the parameter h 
of Eq. (9) for different values of K (from top to bottom: 
K = 0.05,0.1,0.2,0.3,0.4,0.6,1,2,5) and U = 0.6. The re- 
gion in which the parametric resonance occurs is bounded by 
the stars. 



IV. DYNAMICS OF THE FULL SYSTEM 



A. Harmonic case 



Now we consider the original system of equations of 
motion Eq. (5). The linearization given by Eq. (6a) al- 
lows to understand some of the dynamical features of 
the full system. However, the CM and relative motion 
equations are coupled and this results in a qualitatively 
different behaviour with respect to the simple approxi- 
mation discussed in Sect. III. In particular, we note that 
the feedback of x r on xqm drives the CM out of the in- 
stability window found for the linearized system Eq. (9). 
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FIG. 3. Comparison between relative and CM motion of 

the linearized system Eq. (6a) (upper panel) and those ob- 
tained by integrating the complete system Eq. (4) (lower 
panel). The relative motion plotted in (a) is the numerical 
result of the integration of Eq. (9) . The CM and relative co- 
ordinate of the linearized equations, which are plotted in (b), 
oscillate with the same frequency. The CM of the full system 
is shown in (d) where we note a decrease of amplitude at the 
point in which the relative motion starts to decrease in (c). 
Only the envelope of the rapid oscillations of x r is visible on 
the left panel ((a) and (c)). The parameters used in the sim- 
ulation are Uo = 0.6, K = 0.05 and vo = 0.7. All lengths are 
rescaled to the substrate lattice constant a. 
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An example is illustrated in Fig. 3, where we compare 
the motion of the complete and of the linearized system. 
As we can see in Fig. 3(c), the parametric increase of 
x r found for the linearized equations is followed by a de- 
crease, due to the fact that the feedback of x r on xcm 
causes a change of the amplitude of the CM at that point 
(t ~ 400 in Fig. 3(d)). However, when x r decreases, the 
instability reappears and the relative motion increases 
again (although it is not shown in the figure). The sys- 
tem gets in and out the instability window, because we 
are considering a case that is almost at the border. In- 
stead, in Fig. 4 we show a case which is in the centre of 
the region of instability: we can note that after an ini- 
tial transient the relative motion is always excited, and 
its behaviour is more irregular so that it is not possible 
to identify a clear unique frequency of oscillation and a 
unique rate of increase. This behaviour is caused by the 
shift in position inside the instability window which in 
turn produces a shift in frequency and rate of increase. 
Note that the excitation of the internal motion leads to a 
CM motion which would not have occurred if the dimcr 
had been rigid. In that case the CM would have kept 
the initial oscillatory behaviour around the equilibrium 
position. Here, instead, the internal vibrations play a 
role similar to that of a heat bath and drive the CM 
away from the minimum with jumps across one or more 
potential wells. 
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FIG. 4. Numerical simulations of Eq. (5) for Uo = 0.6, 
K = 0.05 and v = 1.54. The CM (solid line) and relative 
motion (dashed line), rescaled to the substrate lattice con- 
stant a, are shown. 

In fact, for a non-rigid dimer, it is possible to get out of 
the well even if vq < \/4Uo- This happens because if the 
internal motion is excited, it is possible that one particle 
remains in the minimum whereas the other reaches the 
nearest maximum. In this way the energy balance is: 

E° kin = \v\ + \v\ = 2U a + \K{a/2f (11) 

and if K is sufficiently small the right-hand side is smaller 
than AUq (we assume v\ = v 2 = v ). Thus vibrational 



energy can be effective in overcoming the barrier. The re- 
sulting motion of the dimer (Fig. 4) can be characterized 
as chaotic, as shown later in Sect. V. 

It is interesting that the chaotic motion described 
above occurs at velocities below and above the thresh- 
old V4J7o for drift motion in the linearized system. This 
is due to a coexistence of long jumps with localized mo- 
tion which persists for a certain range of initial energies, 
as suggested in Ref. [6]. In Fig. 5(a) we show the case 
where v = \/4t/o- At the beginning, the CM performs a 
step-like motion: every time it overcomes a barrier it gets 
stuck for a while in the next minimum before overcom- 
ing the next barrier. In this initial stage x r — 0. After 
x r gets excited, this step-like motion disappears. The 
parametric resonance which one would have expected in 
this case for the simplified system is not visible because 
of the reciprocal influence of CM and relative motion, 
which inhibits the increase of amplitude of x r . Note that 
when the internal motion is excited almost all the energy 
is transferred to the vibrational modes, as it can be seen 
by the corresponding peaks in Fig. 5(d). 




t t 
FIG. 5. Numerical simulations of Eq. (5) for Uo = 0.6, 
K = 0.1 and v = y^Wo ^ 1.5492: (a) CM trajectory (solid 
line) and relative motion (dashed line) ; (b) kinetic energy; (c) 
substrate potential energy; (d) spring potential energy. 

By further increasing the initial kinetic energy the dy- 
namics becomes again non chaotic, and the CM performs 
a drift motion with constant velocity unless the condi- 
tions for parametric excitation given by Eq. (9) are met. 
This does not occur for the small values of K considered 
up to now. In fact, a dimer presents only one charac- 
teristic frequency so that conditions for parametric ex- 
citation are generally met either in the oscillatory or in 
the drift regime. The situation would be different for a 
larger molecule where different vibrational modes could 
be excited for different values of vq. In Fig. 6 we show 
one situation for large K where the drift CM motion ex- 
cites the internal motion through a parametric resonance 
for drift velocity twice the dimer natural stretching fre- 
quency luq — \J2K. When the relative motion acquires a 
large amplitude, deviations from the linear behaviour of 
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xcm are observed. 




FIG. 6. Numerical simulations of Eq. (5) for Uo = 0.6, 
K — 0.5 and v — 2.5. In (a) we show the CM motion (solid 
line) with a linear fit for t < 200 (dashed line), while the 
relative coordinate is plotted in (b). 
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FIG. 7. Numerical simulations of Eq. (5) for Uo = 0.6, 
three values of K and several values of vo- The CM mo- 
tion is plotted for K = 0.05 (a), K = 0.3 (b) and K = 0.5 
(c). The different curves in each graph are obtained with 
different initial velocities. From bottom to top in each 
graph v = 0.63,1.65,1.68 (a), v = 1.5,1.85,2.15,2.2 (b), 
v = 1.54, 1.55, 2.2, 2.5, 2.55 (c). Note the deviation from lin- 
ear behaviour in (c) for vo = 2.5 (see Fig. 6 and text). 

Lastly for K = 0.5 (Fig. 7(c)), there is no chaotic 
regime. The drift motion starts at Vo = y/4Uo and devi- 
ations only occur for a narrow range of velocities where 
x r becomes parametrically excited, as shown in Fig. 6. 

We may estimate from Eq. (11) the critical K value 
above which the internal motion is not effective in making 
the CM overcome the potential barrier for v < ^f4^Q. 
Namely 



In Fig. 7 we summarize for three values of K the ef- 
fect of different vo (at fixed Uo) on the CM motion. In- 
creasing vo, a complex transition from oscillatory regular 
motion to chaotic motion and then to a drift regime can 
take place depending on the value of K. For K = 0.05 
(Fig. 7(a)), resonant excitation of the internal motion 
occurs for v > 0.65 as given in Fig. 1(d). Above this 
value, first a regime with recursive excitation of 
in Fig. 3(c), takes place so that < xcm >= a/2 and 
< vcm >— 0. At u > 1.3 the escape from the well 
described by Eq.( 11) becomes possible. The resulting 
(chaotic) behaviour of xcm in this regime is shown by 
the dashed line in Fig. 7(a) for vq = 1.65. The CM mo- 
tion in this regime behaves as < x 2 CM >~ t a (1 < a < 2) 
with < vcm >— 0, i.e. it is "quasi-diffusive", that is to 
say neither purely diffusive nor ballistic. This behavior 
extends up to v < 1.68 , i.e. well above the "threshold" 
\f£U~o~. Above, a drift motion xcM^t) ~ x + < vcm > t 
occurs. For a larger K = 0.3 (Fig. 7(b)) the quasi- 
diffusive motion starts occurring at values of vq ~ ^/4Uo 
up to uo = 2.2, where drift motion is recovered. 



2U + -K c ir 2 



4U, 



0i 



(12) 



and for Uo = 0.6 we find K c ~ 0.35. As a consequence 
for K > K c no chaotic motion is found. 



B. Lennard- Jones case 

Now we consider the effect of replacing the harmonic 
interaction (2) with a finite-range potential. As a simple 
choice, we take the LJ potential, given by 



V LJ (r) = 4e 



r) 



(13) 



with r = \x2 — X\\ and a cutoff at r — 2.5o\ To recover 
the harmonic interaction close to the minimum, we im- 
pose the equilibrium distance to be equal to the spring 
equilibrium length and the second derivative of Vlj to 
be equal to the spring constant, namely: 



d 2 V L j 



_ 24e 



V2a 

26 _ 7 
vT28 vT6 



= I 

= K 



(14) 
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whence 



{ t 72 



At variance with harmonic interactions, a finite-range po- 
tential allows dissociation of particles. This is clearly seen 
in Fig. 8, where the time behaviour of the CM and rel- 
ative motion is plotted comparing the harmonic and LJ 
potentials. 
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FIG. 8. Numerical simulations of the equation of motion 
of a dimer for harmonic (a),(b) and LJ interactions (c),(d). 
The CM and relative motion are shown. The parameters are 
Uo — 0.6, K = 0.05, vo — 1.12 for the harmonic poten- 
tial. The parameters of the LJ have been chosen according to 
Eq. (15): a ~ 5.598 and e ~ 0.0274. 

The CM is the same in the two cases when x r = 0, i.e., 
at the beginning of the motion. Then, as the relative 
motion starts to increase, x r given by LJ is found very 
similar to the harmonic x r , but then the amplitude of 
the oscillations due to the LJ potential becomes larger. 
At t ~ 420 the two particles dissociate and x r starts to 
increase very fast since only one particle moves. 

In Fig. 9 we show a similar process for the case where 
the CM performs a drift motion. We observe that the 
increase of amplitude of x r occurs approximately at the 
same time for both harmonic and LJ interactions. As just 
noted above, at this point a departure from the linear be- 
haviour of the CM takes place. While large oscillations 
persist in the harmonic x r , breaking of the interparticle 
bond is found in the LJ case. This shows that the res- 
onant excitation of internal vibrations could be effective 
in leading to dissociation of molecular bonds. 
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FIG. 9. Numerical simulations of the equation of motion 
of a dimer for harmonic (a),(b) and LJ interactions (c),(d). 
The CM and relative motion are shown. The dashed lines in 
(a) and (c) are linear fits to xcm for t < 80. The parameters 
are Uo = 0.6, K = 0.4 and vo = 2.15 for the harmonic poten- 
tial. The parameters of the LJ have been chosen according to 
Eq. (15): o ~ 5.598 and e ~ 0.219. 

V. CHAOTIC DYNAMICS 

The dynamics described in Sect. IV A shows very com- 
plex features, in spite of the simplicity of our model. The 
quasi-diffusive, irregular motion found for small values of 
K, as in Fig. 4, resembles characteristics peculiar to a 
chaotic regime. This is confirmed by looking at the tem- 
poral evolution of two trajectories starting at infinitcs- 
imally distant points. For example, Fig. 10 shows two 
long CM trajectories with initial spatial conditions dif- 
fering from 10 -6 . The behaviour of the CM is unpre- 
dictable and the trajectories diverge for the entire sim- 
ulation time. This is a qualitative signature of chaotic 
dynamics. In order to characterize more quantitatively 
the chaotic motion, we have numerically computed the 
Lyapunov exponent, which measures the rate of diver- 
gence of nearby trajectories (see [28,29]): 

8x{t) ~ fa(0)e At , (16) 

where 8x(t) denotes the separation between nearby tra- 
jectories and A is the Lyapunov exponent. 

In an n— dimensional phase space n Lyapunov expo- 
nents can be calculated, but we limit ourselves to the 
computation of the maximal Lyapunov exponent A TOax , 
which is sufficient to signal the occurrence of chaos. If 
A TOO z > the motion is unstable and chaos may oc- 
cur, while if X max < we have a regular stable motion 
(Amarr = corresponds to a stable quasi-periodic mo- 
tion). We show Xmax as a function of time for a small 
value of K (K = 0.05 and U = 0.6) in Fig. 11 for dif- 
ferent vo. The saturation values of the different curves 
give a measure of the corresponding maximal Lyapunov 
exponent. 
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FIG. 10. CM motion for Uo = 0.6, K = 0.1 and 
wo = V^Uo — 1.5492, starting from different initial condi- 
tions: xi(0) = 0, 2:2(0) = a (solid line) and xi(0) = 10 , 
2:2(0) = a + 10~ 6 (dashed line). 

We note that for low values of v (v — 0.63 in the 
figure) Xmax = 0, meaning that the motion in this range 
is regular: xcm oscillates periodically and x r ~ 0. When 
the internal motion starts to be excited (0.65 < v < 1.4) 
^max jumps to a positive small value (X ma x — 0.02), sig- 
nalling that a weak chaotic dynamics is induced by the 
relative coordinate. 
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FIG. 11. Temporal behaviour of the maximal Lyapunov ex- 
ponent for Uo = 0.6, K — 0.05 and different initial velocities 
vo, which are reported beside each curve, except for the two 
curves that saturate at zero, which correspond to vo = 0.63 
and vo = 1.68. 

For larger values of v (1.4 < v < 1.68), as the CM 
gets out of the minimum of the potential well and per- 
forms an irregular motion of the type shown in Fig. 4, the 
magnitude of \ max suddenly increases of about one order 
of magnitude (X ma x — 0.1), but jumps again discontin- 
uously to zero when v is high enough for the CM to 
perform a drift motion and x r ~ 0. In this way, we have 
a complex transition from non-chaotic to chaotic motion 
and viceversa as a function of the initial velocity. This 
behaviour is different from that of the double pendulum, 



where the non-chaotic regime is not recovered for large 
initial velocities [24,25]. 



As a further indicator of such a dynamical behaviour 
we have plotted the phase space projected on the 
(xcm,vcm) plane in Fig. 12, for four different initial 
velocities. 




-5 5 10 15 20 25 50 100 150 200 

X CM X CM 

FIG. 12. Phase space plot projected on the (xcm,vcm) 
plane for Uo = 0.6, K = 0.05 and four values of vo'- vo = 0.63 
(a), v = 0.8 (b), v = 1.4 (c) and v = 1.68 (d). 



The phase plot in (a) is a simple closed loop corre- 
sponding to a regular oscillatory motion where \ max = 0. 
As the initial velocity increases more complex features 
appear: in the weak chaotic regime (b) extra loops are 
present, while the phase plot in (c) becomes very much 
folded and irregular. The regular dynamics is restored 
again in (d), where there is a drift motion of the CM, 
with vcm oscillating around the drift velocity. 



Power spectrum analysis is usually considered as an 
additional effective method to detect chaos. We have cal- 
culated the power spectra of x r by using a fast Fourier 
transform and we show them in Fig. 13 for the same val- 
ues of Uo, K and v as in Fig. 12. The power spectrum 
for the regular motion (a) is smooth and has few peaks, 
at oj osc , SuJosc, 5w osc , i.e. the harmonics expected for 
a parametric oscillator. In (b) each peak broadens, de- 
veloping further lateral features. For the most chaotic 
motion (c) the power spectrum becomes very irregular 
with a large number of peaks. This chaoticity disappears 
for higher velocities (d), where the motion is regular and 
the power spectrum is again smooth with a large peak 
at uj = \J1K corresponding to the dimer stretching fre- 
quency. 
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FIG. 13. Power spectra of the relative coordinate |a; r (ct;)| 2 
for Uo — 0.6, K — 0.05 and four values of vq: vq = 0.63 (a), 
v = 0.8 (b), v = 1.4 (c) and v = 1.68 (d). 

VI. CONCLUSIONS AND DISCUSSION 

In this paper, we have studied a one-dimensional model 
of a dimer moving on a rigid periodic substrate. We have 
shown how the nonlinear effects due to the substrate po- 
tential and to the interplay between the CM and relative 
motion can be relevant in determining the peculiar char- 
acteristics of the dimer motion. A complex dynamical be- 
haviour is found as a function of the initial kinetic energy 
with the occurrence of resonant instabilities and chaotic 
motion. In particular, for weakly bound dimcrs a chaotic 
regime is found for a whole range of velocities between 
two non-chaotic phases at low and high kinetic energy. 
We have characterized this chaotic regime by studying 
the Lyapunov exponents and power spectra. Moreover, 
we have shown that if more realistic, finite-range interac- 
tions are considered, dimer dissociation can be induced 
via this mechanism by choosing appropriate initial veloc- 
ities. 

Although our model neglects thermal fluctuations, we 
can try to make some qualitative predictions concerning 
the effect of a finite substrate temperature, modelled by 
stochastic forces and a damping term, on the diffusive 
dynamics. The introduction of these effects could smear 
out all the deterministic effects. The stable periodic or- 
bits become attracting centres via dissipation, so that the 
regular and chaotic motions would be only transient. On 
the other hand, temperature should be effective to pro- 
vide energy to escape from an attractor, giving rise to a 
diffusive motion. Nevertheless, we expect that for small 
temperatures and small friction coefficients, the thermal 
equilibration time should be bigger than the equilibration 
time due to the deterministic chaotic dynamics [6] . Thus, 
under such circumstances, the effects explained in this 
paper could be significant also for real systems at finite 
temperature at short time scales. In particular, diffusion 
should be highly promoted for weakly bound dimers for 
which we found the chaotic features. Indeed, as reported 
in Ref . [9] , the diffusion coefficient decreases by about an 



order of magnitude with respect to the non-interacting 
case K = 0, when the elastic constant is increased from 
K = to K = 0.25, at least for small values of the 
damping and the temperature. Moreover, signatures of 
the chaotic regime can be recognized in the jump length 
distribution: jumps of the adsorbate over many lattice 
parameters are predicted by theoretical models (see for 
example [11,30]) and observed experimentally (see [31]). 
A preliminary study of the same model at T ^ shows 
that the coexistence of localized and unbounded motion 
of the type shown in Fig. 12(c), is present up to tempera- 
tures a few times higher than the potential barrier 2Uq. 
The relation between chaotic deterministic diffusion and 
stochastic thermal diffusion is an important topic cur- 
rently under study [6,32]. This represents a connection 
between the behaviour of our simplified model and the 
one of more realistic systems. Therefore we believe that, 
beside their intrinsic interest, our results can be of impor- 
tance to understand the dynamical behaviour of dimcrs 
moving onto real surfaces. 

Note added in proof. After completion of this 
manuscript we have noted the recent paper A. S. Ko- 
valev, A. I. Landau, Low Temp. Phys. 28, 423 (2002) 
presenting related numerical results of diffusive dimer dy- 
namics. 
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